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Abstract 

The semi-exponential basis set of radial functions (A.M. Frolov, Physics Letters A 374, 
2361 (2010)) is used for variational computations of bound states in three-electron atomic 
systems. It appears that semi-exponential basis set has a substantially greater poten- 
tial for accurate variational computations of bound states in three-electron atomic systems 
than it was originally anticipated. In particular, the 40-term Larson's wave function im- 
proved with the use of semi-exponential radial basis functions now produces the total energy 
-7.47805413551 a.u. for the ground 1 2 S— state in the °°Li atom (only one spin function xi = 
af3a — (3aa was used in these calculations). This variational energy is very close to the exact 
ground state energy of the °°Li atom and it substantially lower than the total energy obtained 
with the original Larson's 40-term wave function (-7.477944869 a.u.). 
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I. INTRODUCTION 



In this study we perform variational calculations of bound states in three-electron atomic 
systems. An advanced basis set of semi-exponential radial functions [2| is extensively used 
in our calculations. The main goal is to solve the non-relativistic Schrodinger equation 
= E^, where E < and bound state wave function \l/ has the unit norm. The general 
non-relativistic Hamiltonian H of the three-electron atomic problem is (see, e.g., |l|) 



H 
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vi + vi + v 2 3 + -^vj - - — - — —+—+—+— (i) 
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where Ti — is the reduced Planck constant, m e is the electron mass and e is the electric 
charge of electron. In this equation and everywhere below in this study the subscripts 
1, 2, 3 designate the three electrons e~~, while the subscript 4 denotes the heavy nucleus 
with the mass M (M 3> m e ) and positive electric (nuclear) charge Qe. The notations 
rij =| r i — Tj |= Tji stand for the six interparticle distances (= relative coordinates) defined 
in an arbitrary four-body system and r, (i = 1, 2, 3, 4) are the Cartesian coordinates of 
the four point particles. In Eq.flT]) and everywhere below in this work we shall assume 
that (ij) = (ji) = (12), (13), (14), (23), (24), (34). In fact, below only atomic units 
h = 1,| e |= 1, m e = 1 are employed. In these units the explicit form of the Hamiltonian H, 
Eq.([T]), is significantly simplified. 

The main attention in this work is focused on numerical calculations of the ground (dou- 
blet) 1 2 S(L = 0)— state (or l 2 Si(L = 0)— state) of the three-electron Li atom with the 
infinitely heavy nucleus, i.e. the °°Li atom. Considerations of other three-electron atoms, 
ions and various positron containing atomic systems (e.g., HPs) can be performed absolutely 
analogously (for more detail, see [2j) and, therefore, these systems will not be considered 
here. 

The problem of highly accurate calculations of the bound states in three-electron atomic 
systems has attracted a very significant attention. First calculations of the Li atom with the 
truly correlated wave 
tions can be found in 

current bibliography on this subject includes almost one thousand references and rapidly 
increases. At this moment the classical Hylleraas method (see, e.g., [5[) is only the method 
used to construct highly accurate wave functions for bound states in three-electron atomic 



unctions were performed in 1936 jsj]. A brief reviews of such calcula- 
4] (earlier works) and {5] (references up to 1997 are mentioned). The 
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systems. In this method to produce the highly accurate wave functions, e.g., for the ground 
(doublet) l 2 S'-state in the Li atom one needs to use many thousands Hylleraas basis func- 
tions. The use of extremely large basis sets is very inconvenient in many actual cases, since 
it produces a number of computational problems. It is clear that the Hylleraas method can- 
not be used to construct compact and accurate variational wave functions for three-electron 
systems. Indeed, it contains essentially no control parameters which can be optimized by 
increasing the overall efficiency of the method. 

An alternative approach to variational bound state calculations in three-electron atomic 

n 

systems was proposed in |2J. This approach is based on the use of semi-exponential varia- 
tional wave functions ^ and it allows one to construct very compact and accurate variational 
wave functions in arbitrary three-electron atomic system. Each of the semi-exponential ba- 
sis functions depend upon all six interparticle coordinates 7*12 , ris, . . . , Q. A very high 
efficiency of this new approach in actual applications and its superiority over the classi- 
cal Hylleraas expansion was demonstrated in [2]. As follows from the results of this study 
the semi-exponential basis set has a substantially greater potential for highly accurate varia- 
tional computations of bound states in three-electron atomic systems than it was anticipated 
earlier 

II. VARIATIONAL WAVE FUNCTION 

The variational wave function of the doublet S(L = 0)— states of the three-electron Li 
atom is written in the following general form 

*l=o = i>L=a(A; {rij})(a/3a - Pact) + (pL=o{B; {r ij })(2aa/3 - f3aa - a/3a) (2) 

where -0£ =o (^4; {t^-}) and (j)L =0 (B; {r^}) are the two independent spatial parts (= radial 
parts) of the total wave function. Each of these two radial functions is, in fact, a radial 
factor (for states with L = 0) in front of the corresponding three-electron spin functions 
Xx = otfia — (3aa and \2 = 2aa(3 — (5aa — af3a. Here the notations a and f3 are the one- 



electron spin-up and spin-down functions, respectively (see, e.g., js]). The notations A and B 
in Eq. (T5]) mean that the two sets of non-linear parameters associated with radial functions 
and cj) can be optimized independently. In the general case, each of the radial basis functions 
explicitly depends upon all six interparticle (relative) coordinates r 12 , r 13 , r 2 3, r 14 , r 2 4, r 34 . It 
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is clear that in actual variational calculations only one spin function, e.g., the \i function, 
can be used. Another useful trick (so-called 'doubling') is based on the use of the same set 
of non-linear parameters in the two radial parts in Eq.(J2]). 



In our earlier work [2| we have introduced an advanced set of radial basis functions for 
three-electron atomic calculations. In [2| this set was called the semi-exponential basis set. 
In general, the semi-exponential variational expansion of the radial function if)L = o(A; {ry}) 
is written in the form 

N 



ij L=0 (A; { rij }) = £ C fe r 2 "3 l(fc) r- (fc) rr 2 3(fc) rr 4 l(/£) r 2 7 (fc) r 3 7 (fe) e^(-a fc r 14 - /3 fc r 24 - 7fe r 34 ) (3) 

k=l 



where akihilk (k = 1,2, . . . , N) are the varied non-linear parameters. The presence of 
the varied non-linear parameters in Eq.([3]) is the main and very important difference with 
the traditional Hylleraas variational expansion (see, e.g., (7|) for which in Eq.flH]) we always 
have ax = ... = &n, Pi = ■ ■ ■ = 0n and 71 = . . . = 7^. Note that all matrix elements of 
the Hamiltonian, Eq.flT]), and overlap matrix needed in computations with the use of the 
semi-exponential basis, Eq.(j3]), contain the same three-electron integrals which arise for the 
usual Hylleraas expansion (for more detail, see [2f). In other words, numerical calculation of 
all matrix elements with our semi-exponential functions is no more difficult problem, than 
for the traditional Hylleraas radial functions. This also simplifies numerical computation 
of the bound state properties (i.e. expectation values) in the semi-exponential basis set. 
Our algorithms used in calculations of all required matrix elements is based on the Perkins 
formula for three-electron integrals [8] in relative coordinates. Note also that all calculations 
in this work have been performed with the use of standard quadruple precision accuracy (30 
decimal digits per computer word). 

In actual atomic systems any many-electron wave function must be completely an- 
tisymmetric upon all electron variables, i.e. upon all electron spatial and spin vari- 
ables. For three-electron atomic wave function this requirement is written in the form 
^4i23*(l,2,3) = -#(1,2,3), where $ is given by Eq.© and A e is the three-particle (= 
electron) antisymmetrizer A e = e — A2 — A3 — A3 + A23 + A32- Here e is the identity 
permutation, while A? is the permutation of the z-th and j-th particles. Analogously, the 
operator Ajfc is the permutation of the i-th, j-th and k-th particles. In actual computations 
antisymmetrization of the total wave function is reduced to the proper antisymmetrization 



of corresponding matrix elements (for more detail, see, e.g., {2]). Each of these matrix ele- 



merits is written in the form | O \ ^f), where O is an arbitrary spin-independent quantum 
operator which is truly symmetric upon all interparticle permutations. The wave function 
Eq.([2]), contains the two different radial parts ip and <f>. By performing the integration 
over all spin coordinates from here one finds the four spatial projectors V^,V^<j, = and 
presented in j^J. In fact, the explicit form of the and projectors given in j^] 
must be corrected (there is an obvious misprint in the formulas given in [2]]) 

= \ (A 3 - A 3 + A23 - A 32 ) (4) 

'P(/>i> — 2 (-^13 — A3 + A23 — A32) (5) 

For an arbitrary truly symmetric spin-independent operator O each of these four projectors 
produces matrix elements \ O \ ty) of the correct permutation symmetry (for doublet 
states) between all three electrons. The explicit formulas for all matrix elements obtained 
with the radial basis functions, Eq.(j3]), and for three-electron integrals needed in calculations 
can be found in l2j]. 



ILL CALCULATIONS 



Let us apply the semi-exponential variational expansion, Eq.(|3]), to numerical calcula- 
tions of the ground l^-state in the three-electron °°Li atom. In this study we consider the 
two variational wave functions: (a) the wave function which contains 28 radial basis func- 
tions, Eq.(j3]), and (b) the wave function which includes 40 radial basis functions, Eq.(j3]). 
The results (in atomic units) obtained with these two trial wave functions can be found in 
Table I. Tables II and III contain the corresponding radial basis functions, Eq.(|3]), i.e. the 
powers ni(k), ^(fc), n^k), mi(k), rri2{k), m^{k) of six radial variables r^, ri3, r23, r^, r2A, 
and optimized non-linear parameters ctk,Pk,lk- As follows from Table I our variational en- 
ergies obtained for the ground l 2 S'-state in the °°Li atom with the use of semi-exponential 
variational expansion, Eq.([2]), are substantially lower than the corresponding energies deter- 
mined for this state with the same Hylleraas wave function [7[ . Note also that the non-linear 
parameters used in our method (in Eq.(j3])) are constantly varied. Therefore, it is hard to 
say that the total energies obtained in some calculations are 'final'. Formally, based on the 
known convergence rate(s) for our data and by using a few extrapolation procedures we 
can approximately evaluate the limits to which our variational energies will converge, if we 

5 



could perform an infinite number of variations for the non-linear parameters in Eq.(j3]). Such 
limits for the total energies are shown in the fourth column of Table I. These values indicate 
that, e.g., our 40-term variational wave function can produce, in principle, very accurate 
variational energies, if the optimization of non-linear parameters in Eq.flH]) will continue. 

In {2I Larsson proposed a simple (but useful!) trick which allows one to increase the 
overall accuracy of the trial (doublet) wave function. Later this trick was called 'doubling' 
of the wave function and it was used practically in all calculations of the bound doublet 
states in three-electron atomic systems. The idea of doubling is simple and transparent. 
If we already know the radial function constructed for one spin configuration, e.g., for 
Xi = oi(5a — (3aa from Eq.(l2"l), then we can use exactly the same radial basis function for 
another spin configuration \2 — 2ctaf3 — fiaa — aj3a. Formally, it doubles the total number 
of basis functions in the trial wave function. According to the variational principle the 
total variational energy can only decrease during such a procedure. The problem of linear 
dependence of basis functions is avoided in this procedure, since the two spin functions 
Xi = a(5a — (5aa and X2 = 2aa/3 — (3aa — a(3a are independent of each other. In fact, 
for Hylleraas variational expansion the 'doubling' does not work properly, since there are 
obvious linear dependencies between different radial basis functions in those cases when some 
non-linear parameters coincide with each other (for more detail, see [?j]). In semi-exponential 
variational expansion Eq.(j3]) all optimized non-linear parameters are independent of each 
other. Therefore, the coincidence of the pre-exponential factors in Eq.([3]) is not crucial and 
does not mean that such basis functions are linearly dependent. This drastically simplifies 
the actual 'doubling' for Eq.(|3]). The energies obtained with the use of 'doubling' of our 
variational wave functions can be found in the fifth column of Table I. It is clear that the 
non- linear parameters (or parameters in the exponents in Eq.Q) from the second part of 
the total wave function are not optimal, i.e. for all terms which contain basis functions with 
numbers i > 41 the non-linear parameters are not optimal. These 120 (= 40 x3) non-linear 
parameters can be re-optimized and this drastically improves the overall quality of the total 
wave function. For instance, approximate one-step re-optimization of the last 40 non-linear 
parameters in the wave function gives the ground state energy -7.47805456511 a.u., which 
is much better than the 'doubling' energy from Table I. 

As follows from Table I the doubling is not an effective approach for our trial wave 
functions with the carefully optimized non-linear parameters. However, we can modify the 
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original idea of doubling into something new which is substantially more effective in actual 
computations. To illustrate one of such modifications, let us assume that we have con- 
structed 40-term variational wave function, Eq.([n]) which contains 40 x 3 = 120 carefully 
optimized non- linear parameters a%, ji, . . . , a i0 , /3 40 , 740- At the second step of our pro- 
cedure we can add forty additional basis functions with the same pre-exponential factors 



7"23 r"! k \^2^ r Ti^ r ^i^ r Tl^ k \ but slightly different exponents in Eq.Q. In reality, these 



new exponents have been chosen quasi-randomly from three different intervals, e.g., 



where i — 1, 2, . . ., 40 and \\ x // designates the fractional part of the real number x. Small 
deviations of these new exponents from the known 'optimal' values (i.e. from ai,/3i,ji, 
where 1 < i < 40) produce the extended wave function of 'almost optimal' quality. On 
the other hand, even these small differences between exponents allows one to avoid a linear 
dependence between basis vectors in Eq.(j3]). Obviously, this procedure can be repeated a 
number of times. This allows one to construct very accurate trial wave functions which 
contain not only 80, but 400, 800 and even 2000 basis functions with almost 'optimal' non- 
linear parameters. 

It is very interesting to perform variational calculations of the ground state of the °°Li 
atom with the use of 60-term wave function constructed from the analogous 60-term wave 
Larsson's wave function [7|. The variational total energy obtained in (7| with that wave 
function was -7.47801035965 a.u. Our 60-term trial wave function with one spin function 
X, constructed in Q corresponds to the substantially lower total energy -7.478057561 a.u. 
After publication of [2] we have started re-optimization of the non-linear parameters in our 
60-term trial wave function. The current total energy is -7.47805637319 a.u. (only one spin 
function xi is used in our calculations). This energy is above the value from Q]. However, our 
current total energy rapidly decreases with almost constant rate ~ 2.5 • 10~ 7 a.u. per cycle, 
i.e. per one variation of all 180 (= 60 x 3) non-linear parameters in the trial wave function. 
We expect that after an infinite number of variations of the non-linear parameters the total 
energy of the °°Li atom obtained with our 60-term radial function and one spin function 
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will converge to the value -7.4780603(3) a.u. which is very close to the actual ground state 
energy. It will be an outstanding result to obtain the value lower than -7.4780602 a.u. for 
the total energy of the ground l 2 ^— state in the °°Li atom by using only 60-term variational 
wave function. 



IV. CONCLUSION 

The semi-exponential variational expansion [2j is applied for bound state calculations of 
three-electron atomic systems. It is shown that this variational expansion allows one to 
construct very compact and accurate variational wave functions for three-electron atomic 
systems. Currently, the use of semi-exponential radial basis functions is the only way to 
produce compact and accurate wave functions for three-electron atomic systems. The total 
energies obtained in this study for the ground l 2 S'-state of the °°Li atom (see Table I) are 
more accurate than our earlier results from Q and substantially more accurate than the 
original Larsson's wave function [?| with the same number of terms. It indicates clearly 
that our semi-exponential variational expansion, Eq.(j3]), has a substantially greater poten- 
tial for variational bound state calculations in three-electron atomic systems than we have 
anticipated originally j^j. Currently, we continue the process of numerical optimization of 
the non-linear parameters in our trial wave functions constructed with the use of semi- 
exponential variational expansion. Note that the choice of Larsson's wave function(s) as 
the first approximation to the semi-exponential variational expansion is not crucial for our 
method. Many other (different) choices are also possible. 

The semi-exponential variational expansion is also applied for accurate computations of 
various atomic and quasi-atomic three-electron systems. In particular, we have made nu- 
merous improvements in our original computer code j^J. It is clear that the semi-exponential 
basis can be used to construct variational wave functions in the four-electron atomic prob- 
lems. In particular, the semi-exponential variational expansion can be used to obtain very 



compact and accurate wave functions of the ground sing 
of the Be atom and Be-like ions (for more details, see 



et l 1 ^— state and triplet 2 3 S— states 



10] and reference therein). It is 



clear that our method can also be used for rotationally excited states in atomic systems, i.e. 
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for states with L > 1, where L is the electron angular momentum. 
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TABLE I: The total energies E (in atomic units) of the 1 2 S(L = 0)— state in the °°Li atom. N 
designates the number of basis functions used. 

N £(Ref.[4]) £7(Eq.(3)) £ a (Eq.(3)) £(Eq.(3), doubling) 
28 -7.477885105 -7.47803638005 -7.4780368(3) -7.47803657859 
40 -7.477944869 -7.47805413551 -7.4780575(10) -7.47805422814 
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TABLE II: Semi-exponential radial basis functions for N = 28 with optimized exponents. This 
wave function produces the total energy E = -7.47803638005 a.u. for the ground 1 2 S— state of the 
°°Li atom. Only one electron spin-function xi = cxjia — j3aa was used in these calculations. 
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0.340570905705403E+01 
0.182145794896563E+01 
0.276859627151255E+01 
0.286729169596989E+01 
0.278208405003856E+01 
0.373836260739635E+01 
0.337994212648889E+01 
0.167110811817316E+01 
0.30652944727521 1E+01 
0.263305048764256E+01 
0.521025323497647E+01 
0.306893947287392E+01 
0. 186312473874842E+01 
0.454831261352940E+01 
0.295559769377747E+01 
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0.368595304619283E+01 
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0.293743657648302E+01 
0.502952034298549E+01 
0.868992405342042E+01 
0.433237156102395E+01 
0.221163211952812E+01 



P 

0.293295186982955E+01 
0.329011437023099E+01 
0.297627929771076E+01 
0.300017942427107E+01 
0.275168552095737E+01 
0.339167728695617E+01 
0.324225523017272E+01 
0.338626357605054E+01 
0.299446682705661E+01 
0.242503640146648E+01 
0.361243051429800E+01 
0.208144332232383E+01 
0.289718514568496E+01 
0.219828523532609E+01 
0.327296616683230E+01 
0.508008324231496E+01 
0.278359204737422E+01 
0.343630387369804E+01 
0.317945977760753E+01 
0.289825235917765E+01 
0.225117021991187E+01 
0.276113990248958E+01 
0.326870035708320E+01 
0.226578275529031E+01 
0.276772288773810E+01 
0.374899612568756E+01 
0.278286578531691E+01 
0.498938978770980E+01 
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0.771546188231103E+00 
0.316360558089102E+01 
0.668289237018917E+00 
0.637008666805219E+00 
0.645451709771 134E+00 
0.658338285572628E+00 
0.108564332122329E+01 
0.867142603936635E+00 
0.688809932201360E+00 
0.776895925198369E+00 
0.660062464430887E+00 
0.123105191974125E+01 
0.996173732965274E+00 
0.158185067718071E+01 
0.672095727321017E+00 
0.636162427376743E+00 
0.620935860917418E+00 
0.698563354607838E+00 
0.652269908234656E+00 
0.685753232846547E+00 
0.990926897297406E+00 
0.802457922378512E+00 
0.453868977046091E+00 
0.669408905271258E+00 
0.950599883558449E+00 
0.516228416453001E+00 
0.112727477282094E+01 
0.647628325318278E+00 
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TABLE III: Semi-exponential radial basis functions for N = 40 with optimized exponents. This 
wave function produces the total energy E = -7.47805413551 a.u. for the ground 1 2 S— state of the 
°°Li atom. Only one electron spin-function xi = oifia — j3aa was used in these calculations. 
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a 

0.309846257882910E+01 
0.190473561336185E+01 
0.278598083293758E+01 
0.288166643409964E+01 
0.294108078850332E+01 
0.354836463958720E+01 
0.287776340439759E+01 
0.126827836395763E+01 
0.310409040480143E+01 
0.304901894182288E+01 
0.410055124696712E+01 
0.319281393757494E+01 
0.212220234455464E+01 
0.385753710808391E+01 
0.290829017069616E+01 
0.762940571197015E+01 
0.337538643016389E+01 
0.338067126368468E+01 
0.376783902187767E+01 
0.302269266874480E+01 
0.174071943396373E+01 
0.196285150420498E+01 
0.308970260481265E+01 
0.251712242684944E+01 
0.321334220677367E+01 
0.690169324304682E+01 
0.719809975272803E+01 
0.363754290732346E+01 
0.356701987266624E+01 
0.244720102707966E+01 
0.315627784159638E+01 
0.325159760628385E+01 
0.338532643878470E+01 
0.348486218895562E+01 
0.575970395842301E+01 
0.352212817969061E+01 
0.384658211093920E+01 
0.359131298628222E+01 
0.294548846533642E+01 
0.348532012166237E+01 



P_ 

0.291892609218063E+01 
0.366143431469990E+01 
0.282092658543858E+01 
0.284028201233856E+01 
0.281478615927614E+01 
0.298528920333113E+01 
0.309168157398350E+01 
0.388838466426963E+01 
0.297686082083507E+01 
0.222521841071977E+01 
0.320833161427739E+01 
0.245033872966999E+01 
0.270066283806003E+01 
0.320997710808722E+01 
0.343736715561301E+01 
0.406380788113988E+01 
0.290752238041757E+01 
0.336098394507054E+01 
0.307834289691412E+01 
0.371054683906740E+01 
0.316379880955191E+01 
0.280066377711626E+01 
0.294746895743190E+01 
0.408949306371305E+01 
0.364406360008282E+01 
0.367442145267222E+01 
0.307482033524797E+01 
0.256273214565107E+01 
0.340559537814167E+01 
0.271840294504025E+01 
0.302366248489895E+01 
0.330340729374998E+01 
0.298984569514637E+01 
0.311367325997526E+01 
0.393521634949650E+01 
0.366200645034756E+01 
0.412493146244247E+01 
0.426731681125525E+01 
0.329899178470137E+01 
0.454562709500504E+01 



7 

0.845633720338414E+00 
0.157818362413849E+01 
0.679821359567157E+00 
0.651027214862940E+00 
0.660626319649503E+00 
0.672210253693381E+00 
0.732676012007879E+00 
0.696918159081991E+00 
0.712624162776256E+00 
0.707931864133508E+00 
0.685502306391946E+00 
0.131890477152115E+01 
0.102980902194010E+01 
0.149747053647776E+01 
0.656999154522055E+00 
0.643367162248368E+00 
0.757319995998642E+00 
0.689553038540639E+00 
0.692428594575284E+00 
0.259164988133888E+01 
0.135613707594630E+01 
0.135193158043399E+01 
0.815386844474244E+00 
0.114312512643339E+01 
0.854918812691407E+00 
0.515986692127293E+00 
0.136078304480511E+01 
0.168423311171287E+01 
0.824959104823405E+00 
0.895838653318275E+00 
0.903526045564038E+00 
0.897876831570133E+00 
0.936601520312646E+00 
0.975295713010522E+00 
0.781294134452016E+00 
0.347518212792864E+01 
0.976334599279947E+00 
0.125620739151171E+01 
0.769541352183112E+00 
0.894468554048998E+00 
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